Studying Turbulence using Doppler-broadened lines: Velocity Coordinate Spectrum

We discuss a new technique for studying astrophysical turbulence that utilizes the statistics of Doppler-broadened spectral lines. The technique relates the power Velocity Coordinate Spectrum (VCS), i.e. the spectrum of fluctuations measured along the velocity axis in Position-Position-Velocity (PPV) data cubes available from observations, to the underlying power spectra of the velocity/density fluctuations. Unlike the standard spatial spectra, that are function of angular wavenumber, the VCS is a function of the velocity wave number k_v ~ 1/v. We show that absorption affects the VCS to a higher degree for small k_v and obtain the criteria for disregarding the absorption effects for turbulence studies at large k_v. We consider the retrieval of turbulence spectra from observations for high and low spatial resolution observations and find that the VCS allows one to study turbulence even when the emitting turbulent volume is not spatially resolved. This opens interesting prospects for using the technique for extragalactic research. We show that, while thermal broadening interferes with the turbulence studies using the VCS, it is possible to separate thermal and non-thermal contributions. This allows a new way of determining the temperature of the interstellar gas using emission and absorption spectral lines.Comment: 27 page, 3 figures, content extended and presentation reorganized to correspond to the version accepted to Ap

Out-Of-Focus Holography at the Green Bank Telescope

We describe phase-retrieval holography measurements of the 100-m diameter Green Bank Telescope using astronomical sources and an astronomical receiver operating at a wavelength of 7 mm. We use the technique with parameterization of the aperture in terms of Zernike polynomials and employing a large defocus, as described by Nikolic, Hills & Richer (2006). Individual measurements take around 25 minutes and from the resulting beam maps (which have peak signal to noise ratios of 200:1) we show that it is possible to produce low-resolution maps of the wavefront errors with accuracy around a hundredth of a wavelength. Using such measurements over a wide range of elevations, we have calculated a model for the wavefront-errors due to the uncompensated gravitational deformation of the telescope. This model produces a significant improvement at low elevations, where these errors are expected to be the largest; after applying the model, the aperture efficiency is largely independent of elevation. We have also demonstrated that the technique can be used to measure and largely correct for thermal deformations of the antenna, which often exceed the uncompensated gravitational deformations during daytime observing. We conclude that the aberrations induced by gravity and thermal effects are large-scale and the technique used here is particularly suitable for measuring such deformations in large millimetre wave radio telescopes.Comment: 10 pages, 7 figures (accepted by Astronomy & Astrophysics

PyCOOL - a Cosmological Object-Oriented Lattice code written in Python

There are a number of different phenomena in the early universe that have to be studied numerically with lattice simulations. This paper presents a graphics processing unit (GPU) accelerated Python program called PyCOOL that solves the evolution of scalar fields in a lattice with very precise symplectic integrators. The program has been written with the intention to hit a sweet spot of speed, accuracy and user friendliness. This has been achieved by using the Python language with the PyCUDA interface to make a program that is easy to adapt to different scalar field models. In this paper we derive the symplectic dynamics that govern the evolution of the system and then present the implementation of the program in Python and PyCUDA. The functionality of the program is tested in a chaotic inflation preheating model, a single field oscillon case and in a supersymmetric curvaton model which leads to Q-ball production. We have also compared the performance of a consumer graphics card to a professional Tesla compute card in these simulations. We find that the program is not only accurate but also very fast. To further increase the usefulness of the program we have equipped it with numerous post-processing functions that provide useful information about the cosmological model. These include various spectra and statistics of the fields. The program can be additionally used to calculate the generated curvature perturbation. The program is publicly available under GNU General Public License at https://github.com/jtksai/PyCOOL . Some additional information can be found from http://www.physics.utu.fi/tiedostot/theory/particlecosmology/pycool/ .Comment: 23 pages, 12 figures; some typos correcte

Modification of Projected Velocity Power Spectra by Density Inhomogeneities in Compressible Supersonic Turbulence

(Modified) The scaling of velocity fluctuation, dv, as a function of spatial scale L in molecular clouds can be measured from size-linewidth relations, principal component analysis, or line centroid variation. Differing values of the power law index of the scaling relation dv = L^(g3D) in 3D are given by these different methods: the first two give g3D=0.5, while line centroid analysis gives g3D=0. This discrepancy has previously not been fully appreciated, as the variation of projected velocity line centroid fluctuations (dv_{lc} = L^(g2D)) is indeed described, in 2D, by g2D=0.5. However, if projection smoothing is accounted for, this implies that g3D=0. We suggest that a resolution of this discrepancy can be achieved by accounting for the effect of density inhomogeneity on the observed g2D obtained from velocity line centroid analysis. Numerical simulations of compressible turbulence are used to show that the effect of density inhomogeneity statistically reverses the effect of projection smoothing in the case of driven turbulence so that velocity line centroid analysis does indeed predict that g2D=g3D=0.5. Using our numerical results we can restore consistency between line centroid analysis, principal component analysis and size-linewidth relations, and we derive g3D=0.5, corresponding to shock-dominated (Burgers) turbulence. We find that this consistency requires that molecular clouds are continually driven on large scales or are only recently formed.Comment: 28 pages total, 20 figures, accepted for publication in Ap

MYRIAD: A new N-body code for simulations of Star Clusters

We present a new C++ code for collisional N-body simulations of star clusters. The code uses the Hermite fourth-order scheme with block time steps, for advancing the particles in time, while the forces and neighboring particles are computed using the GRAPE-6 board. Special treatment is used for close encounters, binary and multiple sub-systems that either form dynamically or exist in the initial configuration. The structure of the code is modular and allows the appropriate treatment of more physical phenomena, such as stellar and binary evolution, stellar collisions and evolution of close black-hole binaries. Moreover, it can be easily modified so that the part of the code that uses GRAPE-6, could be replaced by another module that uses other accelerating-hardware like the Graphics Processing Units (GPUs). Appropriate choice of the free parameters give a good accuracy and speed for simulations of star clusters up to and beyond core collapse. Simulations of Plummer models consisting of equal-mass stars reached core collapse at t~17 half-mass relaxation times, which compares very well with existing results, while the cumulative relative error in the energy remained below 0.001. Also, comparisons with published results of other codes for the time of core collapse for different initial conditions, show excellent agreement. Simulations of King models with an initial mass-function, similar to those found in the literature, reached core collapse at t~0.17, which is slightly smaller than the expected result from previous works. Finally, the code accuracy becomes comparable and even better than the accuracy of existing codes, when a number of close binary systems is dynamically created in a simulation. This is due to the high accuracy of the method that is used for close binary and multiple sub-systems.Comment: 24 pages, 29 figures, accepted for publication to Astronomy & Astrophysic

Habitable Zones and UV Habitable Zones around Host Stars

Ultraviolet radiation is a double-edged sword to life. If it is too strong, the terrestrial biological systems will be damaged. And if it is too weak, the synthesis of many biochemical compounds can not go along. We try to obtain the continuous ultraviolet habitable zones, and compare the ultraviolet habitable zones with the habitable zones of host stars. Using the boundary ultraviolet radiation of ultraviolet habitable zone, we calculate the ultraviolet habitable zones of host stars with masses from 0.08 to 4.00 \mo. For the host stars with effective temperatures lower than 4,600 K, the ultraviolet habitable zones are closer than the habitable zones. For the host stars with effective temperatures higher than 7,137 K, the ultraviolet habitable zones are farther than the habitable zones. For hot subdwarf as a host star, the distance of the ultraviolet habitable zone is about ten times more than that of the habitable zone, which is not suitable for life existence.Comment: 5 pages, 3 figure

Small Scale Structure at High Redshift: II. Physical Properties of the CIV Absorbing Clouds

Keck HIRES spectra were obtained of the separate images of three gravitationally lensed QSOs (UM 673, Q1104-1804, and Q1422+2309). We studied the velocity and column density differences in CIV doublets in each QSO. Unlike the low ionization gas clouds typical of the interstellar gas in the Galaxy or damped Ly alpha galaxies, the spatial density distribution of CIV absorbing gas clouds turns out to be mostly featureless on scales up to a few hundred parsecs, with column density differences rising to 50 percent or more over separations beyond a few kpc. Similarly, velocity shear becomes detectable only over distances larger than a few hundred pc, rising to 70 km/s at a few kpc. The energy transmitted to the gas is substantially less than in present day star-forming regions, and the gas is less turbulent on a given spatial scale than, e.g., local HII regions. The quiescence of CIV clouds, taken with their probable low density, imply that these objects are not internal to galaxies. The CIV absorbers could be gas expelled recently to large radii and raining back onto its parent galaxy, or pre-enriched gas from an earlier (population III) episode of star formation, falling into the nearest mass concentration. However, while the metals in the gas may have been formed at higher redshifts, the residual turbulence in the clouds and the minimum coherence length measured here imply that the gas was stirred more recently, possibly by star formation events recurring on a timescale on the order of 10-100 Million years (abstract abbreviated).Comment: latex file plus 15 postscript figures (45 pages in total); to be published in the ApJ, June 20, 2001 issu

On the use of fractional Brownian motion simulations to determine the 3D statistical properties of interstellar gas

Based on fractional Brownian motion (fBm) simulations of 3D gas density and velocity fields, we present a study of the statistical properties of spectro-imagery observations (channel maps, integrated emission, and line centroid velocity) in the case of an optically thin medium at various temperatures. The power spectral index gamma_W of the integrated emission is identified with that of the 3D density field (gamma_n) provided the medium's depth is at least of the order of the largest transverse scale in the image, and the power spectrum of the centroid velocity map is found to have the same index gamma_C as that of the velocity field (gamma_v). Further tests with non-fBm density and velocity fields show that this last result holds, and is not modified either by the effects of density-velocity correlations. A comparison is made with the theoretical predictions of Lazarian & Pogosyan (2000).Comment: 28 pages, 14 figures, accepted for publication in ApJ. For preprint with higher-resolution figures, see http://www.cita.utoronto.ca/~mamd/miville_fbm2003.pd

Star Clusters

This review concentrates almost entirely on globular star clusters. It emphasises the increasing realisation that few of the traditional problems of star cluster astronomy can be studied in isolation: the influence of the Galaxy affects dynamical evolution deep in the core, and the spectrum of stellar masses; in turn the evolution of the core determines the highest stellar densities, and the rate of encounters. In this way external tidal effects indirectly influence the formation and evolution of blue stragglers, binary pulsars, X-ray sources, etc. More controversially, the stellar density appears to influence the relative distribution of normal stars. In the opposite sense, the evolution of individual stars governs much of the early dynamics of a globular cluster, and the existence of large numbers of primordial binary stars has changed important details of our picture of the dynamical evolution. New computational tools which will become available in the next few years will help dynamical theorists to address these questions.Comment: 10 pages, 3 figures, Te

Towards absolute scales of radii and masses of open clusters

Aims: In this paper we derive tidal radii and masses of open clusters in the nearest kiloparsecs around the Sun. Methods: For each cluster, the mass is estimated from tidal radii determined from a fitting of three-parametric King's profiles to the observed integrated density distribution. Different samples of members are investigated. Results: For 236 open clusters, all contained in the catalogue ASCC-2.5, we obtain core and tidal radii, as well as tidal masses. The distributions of the core and tidal radii peak at about 1.5 pc and 7 - 10 pc, respectively. A typical relative error of the core radius lies between 15% and 50%, whereas, for the majority of clusters, the tidal radius was determined with a relative accuracy better than 20%. Most of the clusters have tidal masses between 50 and 1000 $m_\odot$, and for about half of the clusters, the masses were obtained with a relative error better than 50%.Comment: 11 pages, 7 figures, accepted for publication in Astronomy & Astrophysic